V 


security  classification  of  thi 


REPORT  DOC 


~1.  REPORT  NUMBER  A  A  I^C  Af\ O  CIPIENT’S  CATALOG  NUMBER 

AFIT/CI/NR  84-62T  M  ‘tUO  _ 

4.  TITLE  (md  Submit)  SE  °F  REPORT  ft  PERIOD  COVEREI 

Timing  Performance  Of  Phase-Locked  Loops  In  THESIS/D/I^W/VHON 

Optical  Pulse  Position  Modulation  Communication  _ 

Systems  «•  performing  org.  report  number 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
ClPIENT'S  catalog  number 


*E  OF  REPORT  ft  PERIOD  COVEREO 


THESIS/O^M^TION 


1 7.  AUTHORfaJ 


Donald  Anthony  Lafaw 


9.  PERFORMING  ORGANIZATION  name  and  address 

AFIT  STUDENT  AT:  University  of  Illinois, 
Urbana-Champaign 

I  I.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

AFIT/NR 
WPAFB  OH  45433 


ft.  CONTRACT  OR  GRANT  NUMBERfeJ 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 


IJ.  NUMBER  OF  PAGES 


14.  MONITORING  AGENCY  NAME  ft  ADDRESS  (It  dlllerml  from  Controlling  Office;  IS.  SECURITY  CLASS,  (ol  (hi*  report; 

UNCLASS 


ISa.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  fof  (hie  Report; 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


ELECTED 
SEP  1  3  1984  1 


[  17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  it  different  from  Report) 


is.  supplementary  notes 


APPROVED  FOR  PUBLIC  RELEASE:  I AW  AFR  190-17 


l^Tn  E.  WOLAVER 
Dean  for  Research  and 
SSarfM  Professional  Development 

v  *  "  f  AC  TT  D  4- *.***>  ***  nrn  , 


AFIT%  Wriqht-Patterson  AFBlOH 


I  19.  KEY  WOROS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 


[20.  ABSTRACT  (Continue  on  reverse  side  H  necessary  and  Identify  by  block  number) 


ATTACHED 


DD  i  jan 1473  EOITION  OF  I  NOV  SS  IS  OBSOLETE 

84  09  13  032- 


UNCLASS 

URITY  CLASSIFICATION  OF  THIS  PAGE  (When  Date  Ent ered; 


TIMING  PERFORMANCE  OF  PHASE-LOCKED  LOOPS 
IN  OPTICAL  PULSE  POSITION  MODULATION 
COMMUNICATION  SYSTEMS 


it 


BY 

DONALD  ANTHONY  LAFAW 
B.S. ,  University  of  Illinois,  1983 


THESIS 

Submitted  in  partial  fulfillment  of  the  requirements 
for  the  degree  of  Master  of  Science  in  Electrical  Engineering 
in  the  Graduate  College  of  the 
University  of  Illinois  at  Urbana-Champaign,  1984 


Urbana,  Illinois 


'f  - 


\ 


ABSTRACT 


— ■  In  an  optical  digital  communication  system,  an  accurate  clock  signal 
must  be  available  at  the  receiver  to  provide  proper  synchronization  with 
the  transmitted  signal.  Phase  synchronization  is  especially  critical  in 
M-ary  pulse  position  modulation  (PPM)  systems  where  the  optimum  decision 
scheme  is  an  energy  detector  which  compares  the  energy  in  each  of  M  time 
slots  to  decide  which  of  M  possible  words  was  sent.  A  timing  error 
causes  energy  spillover  into  adjacent  time  slots  (a  form  of  Intersymbol 
interference)  so  that  only  a  portion  of  the  signal  energy  may  be 
attributed  to  the  correct  time  slot.  This  effect  decreases  the  the 
effective  signal,  increases  the  effective  noise,  and  increases  the 
probability  of  error. 


/  N 

3.  Ip-  this  report,  we  simulate?,  a  timing  subsystem  for  a 

A 

.  satellite-to-satellite  optical  PPM  communication  link.  The  receiver 
employs  direct  photodetection,  preprocessing  of  the  optical  signal,  and 
a  phase-locked  loop  for  timing  synchronization .  The  photodetector 
output  is  modeled  as  a  filtered,  doubly  stochastic  Poisson  shot  noise 
process.  The  variance  of  the  relative  phase  error  is  examined  under 
varying  signal  strength  conditions  as  an  indication  of  loop  performance, 
and  simulation  results  are  compared  to  theoretical  relations.  ^  -  - 
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CHAPTER  1 


INTRODUCTION 


In  an  optical  digital  communication  system,  the  laser  is  modulated 
using  an  appropriate  signaling  format.  In  pulse  position  modulation 
(PPM)  systems,  an  accurate  clock  signal  must  be  available  at  the 
receiver  to  provide  proper  synchronization  with  the  transmitted  signal. 
In  some  systems,  timing  information  is  sent  to  the  receiver  in  the  form 
of  a  special  timing  signal  which  is  superimposed  on  the  transmitted  data 
waveform.  This  timing  signal  is  identified  and  separated  at  the 
receiver  and  used  to  ensure  time  synchronization.  Alternatively,  a 
receiver  could  contain  a  subsystem  to  extract  timing  directly  from  the 
data  stream  and  no  special  timing  waveform  would  be  required,  resulting 
in  a  simpler  and  more  power  efficient  transmitter.  A  synchronous 
subsystem  usually  contains  a  local  oscillator  tuned  near  the  expected 
frequency  of  the  incoming  signal.  Correct  frequency  tuning  of  the  local 
oscillator  is  essential  for  synchronous  reception.  However,  the  phase 
of  the  local  oscillator  must  also  agree  with  that  of  the  incoming 
signal.  Phase  synchronization  is  especially  critical  in  M-ary  PPM 
systems  where  the  optimum  decision  scheme  is  an  energy  detector  which 
compares  the  energy  in  each  of  M  time  slots  in  a  frame  to  decide  which 
of  M  possible  words  was  sent.  A  timing  error  causes  energy  spillover 
Into  adjacent  time  slots  (a  form  of  Intersymbol  interference)  so  that 


only  a  portion  of  the  signal  energy  may  be  attributed  to  the  correct 
time  slot.  This  effect  increases  the  probability  of  error.  The  system 
must  also  have  the  capability  to  track  frequency  displacements  due  to  an 
actual  frequency  difference  between  the  receiver  and  transmitter,  a 
drift  of  the  receiver  oscillator,  or  Doppler  shift,  due  to  motion 
between  the  transmitter  and  receiver,  as  in  satellite  tracking. 

The  effects  of  timing  errors  in  optical  digital  systems  have  been 
extensively  Investigated  [1].  Many  methods  have  been  proposed  and 
employed  for  tracking  systems  including  pulse-edge  tracking,  early-late 
gate  detection,  and  maximum  a  posteriori  (MAP)  estimation.  Phase-locked 
loops  have  also  been  studied  for  use  in  both  optical  on-off  keying  and 
optical  sinusoidal  modulation  schemes  {2],  [3]. 

In  this  report,  we  simulate  a  timing  subsystem  for  a 
satellite-to— satellite  optical  PPM  communication  link.  The  receiver 
employs  direct  photodetection,  preprocessing  of  the  optical  signal,  and 
a  phase-locked  loop  for  timing  synchronization.  The  photodetector 
output  is  modeled  as  a  filtered,  doubly  stochastic  Poisson  shot  noise 
process.  The  variance  of  the  relative  phase  error  is  examined  under 
varying  signal  strength  conditions  as  an  indication  of  loop  performance, 
and  simulation  results  are  compared  to  theoretical  relations. 


Higher-order  PPM  systems  demonstrate  superior  performance  and 
require  less  signal  power  for  a  given  probability  of  error  than  other 
digital  signaling  methods.  However,  use  of  PPM  carries  stringent  timing 
requirements.  Timing  offsets  during  decoding  cause  only  a  portion  of 
the  signal  energy  to  be  included  in  the  signal  count  for  the  correct 
time  slot.  Some  energy  will  spill  over  into  adjacent  slots  causing 
intersymbol  interference.  This  spillover  decreases  the  effective  signal 
power,  increases  the  effective  noise,  and  increases  the  probability  of 
error. 

The  performance  criterion  for  block  encoded  systems  is  the  average 
probability  of  making  a  word  error,  PWE.  From  [4],  the  PWE  for  a  block 
encoded  PPM  signaling  format  with  no  timing  error  is 
exp[-(K  +  MK.  )]  « 

PWE  -  1 - 1 - - - l  Pos(k,Ks  + 

k“l 


k-1 

PosU.l^) 


(2.1) 


M-l 


1(1  +  a)M  -  1] 


where 
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k! 


to  t 


t! 


A 

Pos(k,mv)  £ 


k!  exp(“mv> 


and 

K  ■  average  signal  count 
9 

K.  ■  average  background  count 

D 

M  -  time  slots/word. 


(2.2) 


The  probability  of  bit  error,  PBE,  is  related  to  the  word  error 


probability  by 


PBE 


M 


IM  -  1 


PWE 


(2.3) 


Equation  (2.3)  is  useful  for  comparing  the  performance  of  block  encoded 
PPM  to  other  digital  signaling  techniques. 


2.2  Photodetector  Model 

The  output  current  of  a  photodetector,  iQ(t),  consists  of  a  signal 

current  component,  i  (t),  due  to  the  incident  optical  power  on  the 

s 

detector  and  a  Gaussian  noise  (thermal)  current  component  i  (t),  such 

& 

that 


i„(t)  -  is(t)  +  ig(t).  (2.4) 

The  signal  component  ig(t)  of  (2.4)  is  a  non-stationary ,  inhomogeneous 
shot  noise  process  represented  by  [4] 

N(t) 

i  (t)  -  l  G.  w(t  -  t.)  (2.5) 

8  j-1  J  J 

where 

tj  ■  arrival  time  of  the  jth  photon 
N(t)  -  the  total  number  of  arrival  photons  during  (0,t) 

Gj  -  random,  statistically  independent  detector  gains 
w(t)  ■  response  of  the  detector  to  a  single  electron. 

In  response  to  each  arriving  photon,  the  detector  produces  a  primary 
electron-hole  pair.  The  photon  arrival  times  and  thus  the  electron 


7 


release  times  can  be  modeled  as  a  filtered,  doubly  stochastic  Poisson 
process  with  arrival  rate 

00 

yo  -  +  As  l  p(t "  kTw  ■  ckVM)  (2*6) 


where  the  received  optical  pulse 
square  pulse  of  duration  Tw/M 
signal  counting  rate  defined  by 


shape,  p(t),  is  a  unit-amplitude, 
shown  in  Figure  2.2,  where  ls  is  the 


(2.7) 


and 

PQ  *  optical  power  of  the  received  signal 

nd  **  quantum  efficiency  of  the  detector 

h  «  Planck's  constant 

f  »  the  optical  frequency  of  the  signal 

M  *  the  number  of  time  slots  in  a  word 

ck  “(0»l»2,...,M}is  the  random  data  sequence 

\Q  ■  counting  rate  due  to  dark  current 

T  ■  word  width 
w 

T  /M  ■  time  slot  width, 
w 

Note  that  y t)  is  the  average  rate  of  the  photon  arrival  times.  In  any 
given  time  period  of  T  seconds,  the  probability  that  exactly  N  counts 
are  detected  is  given  by 


P[N,(to,t0+T)] 


A'" 

“nT 


(2.8) 


where 
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(2.9) 


2.3  Signal  Preprocessing 

CorrecC  operation  of  Che  PLL  requires  that  a  dominating  frequency 
component ,  close  to  that  of  the  local  oscillator  in  the  PLL,  be  present 
in  the  input  signal.  This  frequency  component  is  what  the  PLL  "locks” 
onto,  enabling  it  to  precisely  track  the  input  signal.  In  analog  radio 
communication,  for  example,  the  tracked  component  is  the  modulating 
waveform  frequency.  If  we  expect  successful  operation  of  the  PLL  under 
the  proposed  operating  conditions,  we  must  ensure  that  the  PLL  input 
contains  an  adequate  frequency  component. 

First  consider  the  output  from  the  photodetector  described  in 
Section  2.2.  Assuming  a  constant  detector  gain,  a  square  optical  signal 
pulse,  p(t),  and  neglecting  thermal  noise,  the  expectation  of  the 
photodetector  output  conditioned  on  the  random  data  sequence,  c^,  is 

E{iD(t)|ck}  -  w(t)  *  yt)  (2.10) 

where  yt)  is  the  photocount  rate  given  in  (2.6)  and  w(t)  is  the 
detector  impulse  response.  The  unconditioned  expectation  of  the 
detector  output  is 


E <iD(t)  }  -  w(t)  *  <y  t)>  (2.11) 

where 

«o  M-l 

<XR(t)>  •  x0  +  *s  ^  £  p<t  -  kt„  -  JVM)  .  (2.12) 


The  second  term  of  (2.12)  is  the  signal  pulse  shape  averaged  over  the  M 
slots  in  a  word.  Since  Che  received  signal  pulse  shape,  p(t),  is 
square,  the  second  term  of  (2.12)  is  a  constant.  Therefore,  (2.12) 
reduces  to 

\ 

<XR(t)>-Xo+~  (2.13) 

and 

E{iD(t)}  ■  constant  (2.14) 

which  implies 

F[E{ip(t)}]  *  2 *( constant)  6(u>).  (2.15) 

From  this  analysis,  we  conclude  that,  for  a  square  pulse  shape,  there  is 
no  trackable  frequency  component  in  an  expected  detector  output  signal, 
only  a  dc  term  is  present.  Modeling  p(t)  as  a  square  pulse  represents  a 
worst  case  situation.  If  p(t)  has  a  shape  other  than  square,  (2.13) 
would  no  longer  hold  and  other  frequency  components  would  appear  in  the 
output  spectrum.  However,  further  analysis  would  be  required  to 

determine  if  the  frequency  components  generated  would  be  of  sufficient 
amplitude  for  the  PLL  to  track. 

We  now  consider  a  form  of  signal  processing  which  will  result  in  a 
trackable  frequency  component  in  the  spectrum  of  the  signal  before  it  is 
applied  to  the  PLL  input.  Figure  2.3  is  a  block  diagram  of  the  proposed 
detector  and  signal  preprocessor  model.  The  use  of  a  linear  operation 


will  result  only  in  a  convolution  in  the  time  domain  or  multiplication 
in  the  frequency  domain  which  will  change  nothing  but  the  magnitude  of 
the  output  spectrum.  A  nonlinear  operation  is  needed  to  produce  a 
frequency  component  at  the  slot  rate.  A  square-law  device  was  selected 
to  produce  the  nonlinear  effect  because  of  its  simplicity  and  ease  of 
implementation  and  analysis. 

Consider  the  expected  value  of  the  square  of  the  detector  output 
current.  From  Reference  [4]  we  have 

E{iJ(t)}  -  <[h(t)  *  AR(t)]2>  +  h2(t)  *  <XR(t)>  (2.16) 

where  h(t)  is  the  combined  Impulse  response  of  the  detector  and  an 
additional  linear  filter  contained  in  the  preprocessor.  Again,  we 
assume  a  square,  unit-amplitude  signal  pulse,  p(t).  Since  p(t)  is  a 
square  pulse,  <X^(t)>  is  again  a  constant  and  thus  the  second  term  of 
(2.16)  is  also  a  constant.  However,  the  first  term  in  (2.16)  can  be 
written  as 

m  oo 

h(t)  *  X_(t)  -  Xn  /  h(t)  dt  +  X  l  g(t  -  kT  -  c.T  /M)  (2.17) 

K  °  s  k— « 

where 

g(t)  -  h(t)  *  p(t).  (2.18) 

For  computational  convenience  let 
00 

/  h(t)  dt  -  1  .  (2.19) 


In  this  case,  we  obtain 


<[h(t)  *  3^(0]^  -  XQ  +— g- 


+  *s  I  I  <g(t  *  Mw  -  c^T/M) 
J—»  k— • 


•g(t  -  jTw  -  CjTw/M)>  . 


(2.20) 


The  first  two  teras  of  (2.20)  are  constants.  Assuming  the  pulses 
occupying  the  time  slots  are  narrow  and  no  significant  overlap  occurs, 
the  cross  terms  of  the  double  summation  in  (2.20)  are  negligible.  The 
more  slots  per  word,  the  better  this  approximation  will  be.  Use  of  this 
approximation  implies  that  the  only  significant  terms  in  the  summation 
are  those  for  which  j»k,  and  we  obtain 


\i  i  <•>•»£! 

j  *—oo 


(2.21) 


so  that 


l  <8Z(t  ‘  kTw  '  clcTw/M)>  *  Xl  l  i  g2(t  '  nT„/M)  *  (2'22) 

k»-«®  n— » 


After  substituting  (2.22),  (2.21),  and  (2.20)  into  (2.16),  we  obtain 


tjco)  -  L  +^.l/'h2(t)  dc  +  1  h(t)  dt 


14 


Since  all  but  the  last  term  are  constants,  (2.23)  is  equivalently 
expressed  as 


E{iJ(t)}  *  K  +—■  [h(t)  *  p(t)]: 


*  l 

k— 


«(t  -  kTw/M) 


(2.24) 


where  K  denotes  the  sum  of  the  constant  terms  of  (2.23). 

Taking  the  Fourier  transform  of  (2.24)  yields 

2  2"*2, 

[EUi(t)}]  -  2*K«(u>)  +-~  [H(  w)P(  u>)  ]  *  [H(u)P(u)] 

D  w  (2.25) 

CO 

x  y  6(u  -  k2nM/T  ) 
k-—  w 


where  H(ui)  and  P(w)  are  the  Fourier  transforms  of  h(t)  and  p(t), 
respectively.  P(u)  is  simply  the  transform  of  the  transmitted  square 
puls'*  and  is  given  by 


T  bfl* 

v  L  M  (JT./2M 


(2.26) 


From  (2.25)  it  is  apparent  that  we  have  spectral  frequency  components  at 


2  nM 

0)  m  — —  n 
n  T 

w 


(2.27) 


When  n-1,  corresponds  to  the  time  slot  transmission  rate.  H(<i>) 
should  be  chosen  to  maximize  the  signal  component  at  and  to  minimize 
the  other  background  noise  components. 

The  next  step  is  to  put  a  constraint  on  the  choice  of  the  combined 
receiver  filter  impulse  response.  First,  the  detector  output  can  be 
represented  by 


i(t)  -  i(t)  +  &i(t) 


(2.28) 


where  the  subscripts  are  dropped  for  convenience,  and  i(t)  ”  E{i(t)}. 
The  squared  output  of  the  detector  is 

i2(t)  -  TTtT2  +  2lTt T  M(t)  +  Ai2(t)  .  (2.29) 

Taking  expectations  on  both  sides  of  (2.29)  and  assuming  E{Ai(t)}  »  0 
give 

i2(t)  -  UtT2  +  M2(t)  .  (2.30) 

Substituting  (2.30)  into  (2.29)  yields 

i2(t)  -  i2(t)  +  2i(tT  Ai( t )  +  I di2(t)  -  Ai2( t) ]  .  (2.31) 

In  comparision  to  the  first  two  terms  of  (2.31),  the  third  term, 
2  o 

&i  (t)  -  Ai  (t)  ,  is  considered  small  and  can  be  neglected  for  a  first 
approximation,  leaving 

i2(t)  *  i2( t)  +  2lTt7  Ai(t)  .  (2.32) 

To  examine  the  power  density  spectrum  of  (2.32),  we  first  calculate 
the  autocorrelation  function 

<i2(t.),i2(t9)>  -  i2(t.)  i2(t,)  +  4i(t.)  i"(t,y 

1  1  1  _  _I _  (2.33) 

x  IKtj)  -  i(t1)Jli(t2)  -  i(t2)]  . 


The  first  term  of  (2.33)  is  the  expected  output  signal  from  the  squaring 
loop  which  contains  a  nonzero  frequency  component;  the  second  term  is 
noise.  To  minimize  the  noise  term,  remember 


model 
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CHAPTER  3 


DESIGN  AND  ANALYSIS  OF  THE  PHASE-LOCKED  LOOP 

3.1  Mathematical  Model  of  the  Phase-Locked  Loop 

A  block  diagram  of  a  phase-locked  loop  (PLL)  is  given  in  Figure  3.1. 
In  the  basic  operation  of  the  PLL,  a  phase  detector  measures  the  phase 
difference  between  its  two  inputs  and  then  outputs  a  voltage 
proportional  to  this  difference,  vd(t).  The  difference  voltage  is  then 
passed  through  a  loop  filter  and  fed  into  a  voltage  controlled 
oscillator  (VCO)  which  changes  its  frequency  in  a  direction  to  minimize 
the  phase  difference  between  the  VCO  output  and  the  input  signal.  The 
loop  is  said  to  be  ''locked"  when  the  control  voltage  applied  to  the  VCO 
keeps  the  frequency  of  the  VCO  equal  to  the  average  frequency  of  the 
incoming  signal. 

For  the  analysis  of  the  PLL,  assume  that  the  input  signal  has  phase 

0  (t),  the  VCO  output  has  phase  0  (t),  and  that  the  loop  is  locked.  The 
s  o 

phase  detector  is  modeled  as  a  perfect  multiplier  of  gain  K^  whose 

output,  over  a  specified  operating  range,  is  a  direct  voltage,  vd(t), 

that  is  a  function  of  the  phase  difference  0  (t)  -  0  (t).  Under  locked 

s  o 

conditions  the  phase  difference  will  be  small  and  the  phase  detector 
output  is  [5] 


where 


K  A  A 
m  o  s 


(3.2) 


is  the  phase  detector  gain  and  Aq  and  Ag  are  the  amplitudes  of  the  VCO 
output  and  input  signals,  respectively. 

To  reject  high  frequency  signal  components,  the  phase  detector 
voltage  is  filtered  by  a  low-pass  loop  filter  of  transfer  function  F(s). 
The  resulting  control  voltage  is  then  applied  to  the  VCO. 

The  VCO  has  a  free-running  frequency  and  an  instantaneous 
frequency  shift  proportional  to  the  control  voltage  input,  A(tf*KQvc(t), 
where  Kq  is  the  VCO  gain  constant  and  v£(t)  is  the  filter  output 
voltage.  The  VCO  output  frequency  is  written  as 


“o  *  “f  +  *  “f  +  K0vc(t)  • 


(3.3) 


Noting  that  frequency  is  the  derivative  of  phase,  we  see  that 


d0o(t) 

dt 


Kovc(t)  * 


(3.4) 


Using  the  Laplace  tranforms  of  (3.1)  and  (3.4)  we  arrive  at  the 
linearized  mathematical  PLL  model  of  Figure  3.2.  The  loop  transfer 
function  of  Figure  3.2  is 


Vc>  ,  W!s> 

8  <t)  "  H(9)  "  a  +  K  I.Ki)  ' 
s  o  a 


(3.5) 


For  a  further  analysis  we  need  to  specify  F(s).  Using  a  simple 


low-pass  filter,  F(s)  is  of  the  form 
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F(s) 


1 

ST  +  1 


(3.6) 


where  1/t  is  the  3  db  cutoff  frequency  of  the  filter.  Substituting 
(3.6)  into  (3.5),  the  loop  transfer  function  becomes 


H(s) 


Vd/T 

s2  +  s/t  +  K  K./t 
o  a 


(3.7) 


From  control  theory  let  us  define  £  as  the  damping  ratio  and  as  the 
natural  undamped  frequency.  The  denominator  of  (3.7)  is  the 
characteristic  equation  of  the  system  that  determines  the  system 
response  and  is  defined  by 

A  -  s2  +  2£<^s  +  (3.8) 


which  implies 


KoKd 


Vi 


(3.9) 


and 
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(3.10) 


The  selection  of  the  damping  ratio  is  of  particular  importance.  With 
too  large  a  value  of  £,  the  system  response  is  overdamped,  resulting  in 
a  very  long  response  time  to  an  input.  With  too  low  a  value,  the  system 
is  underdamped,  causing  oscillation  about  the  steady-state  value  and 
requiring  a  long  time  for  a  settled  final  response. 

Two  final  design  parameters  to  consider  are  the  lock-in  range  and 


1 


2 


hold-in  range,  Ao^  and  Ao^,  respectively.  When  the  difference  between 
the  free-running  VCO  frequency  and  the  average  incoming  frequency  is 
less  than  the  3  db  loop  bandwidth,  the  loop  will  lock  up  almost 
immediately  without  skipping  cycles.  The  maximum  frequency  for  which 
this  is  possible  is  the  lock-in  range  and  is  roughly  equivalent  to  the 
loop  bandwidth  [3].  The  range  of  frequency  difference  over  which  a  loop 
will  maintain  phase  lock  is  the  hold-in  range  and  is  defined  by  [5] 

Ac^  -  !KoKdF(s)  .  (3.11) 

More  detailed  discussion  and  derivation  of  PLL  design  parameters  can  be 
found  in  [5j  and  [6]. 

3.2  Selection  of  Design  Parameters 

Successful  PLL  design  represents  a  compromise  between  performance 
features.  A  wide  loop  bandwidth  will  result  in  a  fast  acquisition  and 
lock-up  time  as  well  as  an  extended  lock-in  range,  but  will  admit  more 
noise,  may  have  a  smaller  hold-in  range,  and  be  more  sensitive  to  small 
changes  in  the  input  signal  frequency.  Alternatively,  a  narrow  loop 
bandwidth  insures  good  noise  performance  and  will  tolerate  a  larger 
amount  of  input  frequency  change  without  losing  lock,  but  will  have 
difficulty  initially  acquiring  phase-lock.  Advanced  techniques  such  as 
frequency  sweeping  and  loop  bandwidth  changing  have  been  developed  for 
acquiring  lock  even  with  very  narrow  bandwidths,  but  they  are  beyond  the 
scope  of  this  report.  A  thorough  treatment  may  be  found  in  [6]. 


The  simulations  described  in  this  report  used  one-second  time  slots 

which  correspond  to  an  expected  input  frequency  (slot  rate)  of 

2ir  rad/sec.  The  VCO  free-running  frequency  was  selected  as 

6.18  rad/sec.  The  difference  of  0.103  rad/sec  represents  a  1.61 

deviation  from  the  expected  input  frequency,  and  is  used  to  account  for 

uncertainties  about  the  exact  frequency  of  the  received  signal.  Such 

uncertainties  are  commonly  caused  by  Doppler  shifts,  drift  in  the 

transmitted  clock,  or  drift  of  the  receiver  local  oscillator.  The  phase 

detector  gain,  VCO  gain  constant,  loop  filter  bandwidth,  and  damping 

ratio  were  chosen  as: 

Kd  -  0.3333  V/rad 

K  *  1.0  rad/sec/V 

o 

1/t  -  0.3333  rad/sec 
C  *  0.5 

which  yield  the  parameters: 

Any  -  0.3333  rad/sec 
Au^  -  0.42  rad/sec. 

It  is  Important  to  recall  from  (2.2)  that  the  phase  detector  gain  Kj 

is  dependent  on  the  amplitude  of  the  input  signal,  A  .  In  fact,  we  are 

s 

implicitly  assuming  that  the  average  magnitude  of  Ag  is  equal  to  one. 
If  the  value  of  changes  due  to  a  change  in  Ag,  the  design  parameters 
of  the  PLL  also  change.  For  example,  if  Ag  increases  from  1  to  10,  the 
3  db  loop  bandwidth  increases/  from  0.42  to  1.6  rad/sec  causing  the  loop 
to  reject  less  noise  than  originally  intended.  Realistically,  the  input 
signal  will  vary  over  a  substantial  range  of  values  and  the  square-law 


device  in  Che  preprocessor  compounds  Che  problem.  To  compensaCe  for 
varying  signal  amplicudes,  an  accenuacor  is  placed  between  che  squaring 
loop  and  Che  PLL  inpuC.  This  device  will  adjusC  Che  preprocessor  ouCpuC 
signal  co  insure  Che  average  PLL  inpuC  signal  will  not  exceed  unicy. 

To  choose  Che  accenuacion  constant,  we  first  must  find  Che  expected 
value  of  Che  amplitude  out  of  Che  preprocessing  components.  As  in 
Chapter  2,  let  g(t)  denote  the  convolution  of  Che  input  pulse  shape, 
p(t),  and  the  combined  receiver  impulse  response,  h(t).  Let  each  puls» 
have  amplitude  PQ  and  width  T*Tw/M  sec.  From  (2.33),  the  receiver 
Impulse  response  is 

h(t)  -  sin  t,  ,  .  (3.12) 


Convolving  h(t)  with  p(t)  yields 
p  t 

-  1  -  cos  t]  ,  -T  <_  t  ±  0 

g(t)  -  (3.13) 

p  t  (  > 

+  1  -  cos  t  j  ,  0  <_  t  _<  T 


Figure  3.3a  shows  that  as  a  result  of  the  convolution  operation,  g(t) 
has  duration  2T.  After  passing  g(t)  through  the  squaring  loop,  we 


obtain 


3pV  pV  ,  P*T2 

g2(t)  .  — 2_ - 2^  co.,  lit  +-2r 

8ir  2  ir  1  8  it 


4  7T 

cos  y~  t. 


-T  <  t  <  T 


(3.14) 


which  consists  of  two  Identical,  positive  waveforms  each  of  duration  T, 

as  shown  in  Figure  3.3b.  Thus,  the  total  amplitude  of  the  fundamental 

2 

frequency  component  of  g  (t)  is 


Figure  3.3b.  Result  of  passing  g(t)  of  Figure  3.3a  through  the  squaring 
loop. 
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(3.15) 


Averaging  (3.15)  over  Che  M  time  slots  In  a  word  and  substituting 
T  *  T^M  give  us  an  amplitude  coefficient  of 


2  2 
PI 
o  w 

M3w2 


(3.16) 
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CHAPTER  4 

COMMUNICATION  SYSTEM  LINK 


Of  importance  in  any  practical  communication  link  is  the  amount  of 
signal  power  received.  The  larger  the  received  signal,  the  better  we 
can  determine  the  original  message  and  the  lower  the  probability  of 
making  an  error.  In  this  analysis,  we  are  ultimately  interested  in  the 
number  of  photocounts  per  signal  pulse  or  time  slot  received  by  the 
photodetector.  Specifically,  we  now  consider  optical  transmission 
between  two  satellites,  one  in  a  geosynchronous  orbit  and  the  other  in  a 
lower  earth  orbit. 


4.1  Transmitter  Design 

The  purpose  of  an  optical  transmitter  is  to  focus  as  much  of  the 
modulated  signal  power  (led  or  laser)  as  possible  onto  the  receiving 
telescope.  For  a  laser  transmission  system,  the  transmitter  usually 
consists  of  a  beam  expander  which  is  used  to  reduce  the  beam  divergence. 
The  beam  expander  may  consist  of  reflecting  or  refracting  optics 
depending  on  cost,  weight,  and  wavelength  of  the  source.  Figure  4.1 
shows  two  beam  expander  configurations. 

The  output  of  the  beam  expander  can  be  considered  as  a  uniformly 
illuminated  circular  diffracting  aperture  of  diameter  d^.  The  amplitude 
distribution  of  the  Fraunhofer  (far-field)  diffraction  pattern  is  17) 


CASSEGRAINIAN  TELESCOPE 


Figure  4.1.  Typical  beam  expander  configurations 


(4.1) 
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exp(jkz)  exp( jkr^/2z) 


2J1(kdTr0/2z) 


kdrro/2z 


where 

rQ  »  radial  coordinate  in  the  far-field  plane 
d^  ■  diameter  of  the  transmitting  aperture 
k  *  wave  number  of  the  source 
z  *  transmission  distance 
and  Jx(x)  is  a  first-order  Bessel  function, 
intensity  distribution  is 


Kro)  -  U(rQ)  U*(rQ) 


2' 
_ T 

'2J1(kdTro/2z) 

hz 

kdTr£>/2z 

Thus 


the  far-field 


(A. 2) 


where  IT  is  the  transmitted  intensity.  This  intensity  distribution  is 
commonly  known  as  the  Airy  pattern  and  has  well-defined  zeros  wherever 
the  quantity 

[2J1(kdTr0/2z)| 


kdTro/2z 


(4.3) 


For  large  distances,  the  receiver  aperture  is  small  compared  to  the 
size  of  the  Airy  disk  (the  circular  area  contained  within  the  first 
minimum  of  (4.3)).  Therefore,  we  make  the  approximation  for  an  aperture 
of  diameter  dr,  centered  at  a  point  rQ  in  the  receiving  plane,  that  the 
total  optical  power,  P^,  intercepted  by  the  aperture  is 


wd 

r 

~T~ 


Ar  I(rQ) 


(4.4) 


where  Af  ■  irdr/4  is  the  receiver  aperture  area.  Also,  assuming  uniform 


Illumination  of  the  transmitting  aperture,  the  transmitter  laser  power 
PT  is  related  to  the  transmitted  intensity  by 


At1! 


where  A^,  *  itd^/A  is  the  transmitter  aperture  area. 
(4.4),  and  (4.5),  we  write 

'2 


Vr 

2Jl(wdTr0/zX) 

z2  X2 

ndj,rQ/z\ 

(4.5) 

Now  using  (4.2), 


(4.6) 


where  we  have  used  k*2ir/X,  where  X  is  the  wavelength  of  the  source. 

We  now  consider  optical  power  losses  in  the  transmitter.  First,  no 
optical  surface  will  transmit  all  of  the  incident  optical  power. 
Compounding  this  problem  is  the  fact  that  there  is  generally  more  than 
one  lens  or  mirror,  and  the  combined  effects  are  multiplicative. 
Second,  since  the  laser  cross-section  can  be  mathematically  described  as 
a  Gaussian  beam  which  is  Infinite  in  extent,  we  must  account  for 
truncation  of  the  expanded  beam  by  the  finite  diameter  transmission 
aperture.  Finally,  if  a  Cassegrainian  or  similar  type  of  reflecting 
telescope  is  used,  the  secondary  mirror  obstructs  the  beam  and  reduces 
the  overall  transmissivity  of  the  telescope.  In  some  cases,  the 
secondary  mirror  mounting  supports  may  also  contribute  to  obscuration  of 
the  beam.  A  detailed  analysis  of  beam  truncation  and  obscuration 
effects  can  be  found  in  [8]  and  [9]. 

For  this  analysis,  we  account  for  all  the  transmitter  losses  by 
defining  a  constant,  r^,  which  is  the  overall  efficiency  (transmittance) 
of  the  transmitting  telescope.  Thus,  (4.6)  becomes 


4.2  Receiver  Design 

The  objective  of  the  optical  receiver  is  to  focus  the  light  gathered 
by  the  aperture  onto  a  photodetector  while  rejecting  as  much  background 
radiation  as  possible.  A  typical  receiver  is  shown  in  Figure  4.2.  The 
field  stop  iris  is  used  to  narrow  the  field  of  view  of  the  receiver  so 
radiation  from  areas  adjacent  to  the  transmitter  is  ignored.  The 
Interference  filter  eliminates  wavelengths  which  are  not  the  same  as 
that  of  the  source.  In  addition  to  the  optical  losses  mentioned  in 
Section  4.1,  the  receiver  has  the  additional  loss  associated  with  the 
interference  filter.  Narrow  band  filters  have  small  transmittances 
whereas  wider  band  filters  have  higher  transmittances.  For  our 
analysis,  as  in  the  transmitter,  we  define  a  constant,  r^,  to  be  the 
overall  receiver  efficiency  to  account  for  all  losses. 

The  receiver  focuses  the  optical  power  onto  a  photodetector  such  as 
an  avalanche  or  PIN  photodiode.  In  turn,  the  incident  photons  cause  the 
generation  of  electron-hole  pairs  which  flow  to  produce  current  in  the 
device.  The  detector  quantum  efficiency,  nd,  is  defined  as 

Number  of  electron  hole  pairs  generated 
^d  ”  Number  of  incident  photons  *  (4*8) 

A  practical  photodiode  may  have  a  quantum  efficiency  of  between  30%  and 
90%.  However,  a  high  quantum  efficiency  requires  a  thick  depletion 


region  of  material  which  results  in  a  slow  response  time 


A  fast 
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response  photodiode,  necessary  for  high  data  rate  communication  systems, 
will  typically  have  a  much  lower  quantum  efficiency  due  to  the  tradeoff 
between  response  time  and  efficiency. 

We  are  interested  in  determining  the  number  of  photocounts  per 
signal  pulse  generated  by  the  photodiode.  Letting  D  denote  the  data 
rate  of  the  signal,  we  write 


<N> 


nd\X 

Dhc 


P 

r 


(4.9) 


where 

<N>  •  expected  number  of  photocounts  per  signal  pulse 

h  ■  Planck's  constant 

c  «  speed  of  light. 

Combining  (4.9)  and  (4.7)  gives  us 

)2 

<N> 


"d'VV^r 

2J1(ndTrQ/zX) 

Dhcz2  X 

mlj^/zX 

(4.10) 


Equation  (4.10)  is  the  total  link  equation  for  the  system.  Note  that 
(4.10)  is  maximized  if  rQ  ■  0  which  corresponds  to  zero  pointing  error 


between  the  transmitter  and  receiver.  So 
ndnrnTATAr  „ 


<N> 


max 


Dhcz2 X 


(4.11) 


is  the  maximum  number  of  photocounts  per  signal  pulse. 


4.3  Background  Radiation 

For  two  satellites  in  an  earth  orbit,  background  radiation  (noise) 


* 


sources  include  the  stars,  planets,  earth,  moon,  and  sun.  Reception 

with  the  sun  directly  in  the  background  is  usually  not  attainable 

because  of  its  high  brightness,  and  will  not  be  considered  here.  A 

useful  radiation  measurement  is  the  spectral  radiant  emittance,  W(X), 

which  is  the  radiant  power  at  a  wavelength  X  emitted  into  a  hemisphere 

per  unit  area  of  source  in  the  hemisphere.  Assuming  a  Lambertian  source 

of  diameter  d  at  distance  z  from  the  receiver,  the  total  amount  of 
s 

background  radiation  power  incident  on  the  detector,  P^,  is 

,  j2  ,2 

d_d_ 

p  _jLis_r  w(x)  (4.12) 

B  16z2 

where 

nr  *  receiver  efficiency 

dg  -  source  aperture  diameter  (m) 

df  -  receiver  aperture  diameter  (m) 

z  *  distance  between  source  and  receiver  (m) 

*  interference  filter  bandwidth  (A). 

Equation  (4.12)  assumes  that  the  whole  source  is  contained  in  the 

receiver  field  of  view.  If  the  solid  angle  of  the  receiver  field  of 

view,  0r,  is  less  than  the  solid  angle  subtended  by  the  source  at  the 

receiver,  8  ,  the  background  radiation  power  is  given  to  a  good 
s 

approximation  by  [10] 
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(4.13) 


The  solid  angle  for  the  receiver  (assuming  small  angles)  is 


which  Is  the  formula  for  an  extended  source  filling  the  receiver  field 
of  view,  as  is  the  case  when  a  satellite  is  pointed  towards  the  earth 
for  signal  reception. 

For  smaller  sources  such  as  stars,  it  is  usually  more  convenient  to 
make  background  radiation  measurements  in  terms  of  the  spectral 
irradiance,  H(X),  which  is  the  power  per  unit  wavelength  interval 
incident  on  a  unit  area  of  the  receiver  [10].  Therefore,  if  a 


background  radiation  measurement  is  given  in  terms  of  H(  X) ,  then  the 


background  radiation  power  gathered  by  the  receiver  is 


'Vidr 


H(  X)  . 


(4.17) 


Of  primary  concern  is  the  background  radiation  from  the  earth  due  to 
reflection  of  the  sun.  Figure  4.3  shows  the  spectral  radiant  emittance 
of  the  earth  versus  wavelength.  Note  that  for  the  wavelengths  of 
interest  (between  0.8  and  0.9 un) ,  W( X)  is  near  its  peak  value.  As 


Indicated  in  Figure  4.3,  this  measurement  is  valid  when  there  is  no 
cloud  cover  in  the  earth’s  atmosphere.  Radiance  from  sunlit  clouds  is 


Spectral  radiant  emittance,  (M. watts  per  cm2,  micro 


approximately  one  order  of  magnitude  larger  and,  as  we  will  see  in 
calculations  to  follow,  it  represents  the  strongest,  worst  case 
background  noise  source.  Figure  4.4  shows  the  spectral  irradiance  for 
the  planets  and  moon,  and  Figure  4.5  shows  the  spectral  Irradiance  for 
the  brighter  stars  immediately  outside  the  earth's  atmosphere. 
Figures  4.3  and  4.4  show  that  a  full  moon  exibits  the  largest  background 
emission  next  to  that  of  the  earth.  Because  of  the  close  proximity  of 
the  satellite  to  the  moon,  the  use  of  (4.17)  would  be  incorrect  since  it 
assumes  the  whole  source  is  in  the  field  of  view,  when  actually  only  a 
portion  of  the  moon  would  be  in  the  satellite's  field  of  view.  Assuming 
that  the  moon  appears  as  a  uniformly  illuminated  Lambertian  surface,  we 
have  the  relation  [10] 

HU)  .  (4.18) 


Substituting  (4.18)  and  (4.15)  into  (4.16)  yields 


"ViW 


(4.19) 


which  is  the  background  power  gathered  by  the  receiver  when  the  source 
fills  the  field  of  view  and  the  spectral  irradiance  is  the  given 
measurement  parameter.  We  can  convert  the  background  radiation  power 
gathered  by  the  receiver  aperture  into  background  counts/ time  slot 
output  by  the  detector  using 


<n>b 


nd* 

DhcM  PB 


(4.20) 


where 


atmosphere.  (From  W.  K.  Pratt,  Laser  Communication  Systems.  New  York 
John  Wiley  and  Sons,  1968,  p.  123.) 


<N>g  -  expected  number  of  background  counts/time  slot 

and 

H  *  number  of  time  slots/word. 

4. A  Calculation  of  Signal  and  Background  Counts 

In  Sections  4.2  and  4.3  we  derived  the  mathematical  relationships 
governing  the  amount  of  signal  and  noise  counts  emitted  by  the 
photodetector  given  the  signal  (laser)  and  background  power  present. 
Now  we  use  some  typical  component  parameters  to  calculate  numerical 
results. 

The  following  parameter  values  were  chosen  for  a  typical 
satellite-to-satellite  link: 


transmission  wavelength 

X  *  830  no 

transmitter  efficiency 

hj  *  0.7 

receiver  efficiency 

“  0.6 

receiver  planar  field  of  view 

8  ■  1  mrad 
r 

receiver  aperture  diameter 

df  ■  0.25  m 

transmitter  aperture  diameter 

dp  ■  0.1  m 

laser  power 

■  200  mW 

quantum  efficiency 

nd  -  10% 

data  rate 

D  ■  2  X  10®  words/sec 

separation  length 

z  -  3.57  X  107  m 

interference  filter  bandwidth 

-  20  A. 

Using  the  above  parameter  values  in  (4.11)  yields 


77  counts/signal  pulse 
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This  Is  the  maximum  number  of  signal  counts  we  can  expect  under  Ideal 
(no  pointing  error)  conditions.  Note  that  pointing  error  tolerances  are 
severe  since  the  Airy  intensity  pattern  falls  away  rapidly  from  the 
central  maximum.  In  fact,  a  pointing  error  of  only  S  urad  causes  <N>  to 
fall  to  29  counts/signal  pulse. 

Using  (4.16)  and  (4.20)  to  calculate  background  counts  due  to  solar 
radiation  reflected  from  the  earth  filling  the  receiver  field  of  view 
gives 

<N>g  *  3.08/M  counts/time  slot  (no  clouds). 

As  previously  discussed,  if  extensive  cloud  cover  were  present  on  earth, 
this  value  would  increase  to  approximately  30.8/M  counts/time  slot. 

The  moon  is  the  other  major  source  of  background  radiation. 
Figure  4.4  shows  that  the  spectral  irradiance  of  a  full  moon  is  roughly 
10  ^  W/m^/um.  Substituting  this  value  into  (4.19)  and  (4.20),  the 
maximum  number  of  background  counts  possible  from  a  full  moon  is 

<N>g  »  1.52/M  counts/time  slot. 

This  number  of  counts  is  attained  only  if  the  receiver  is  pointing 
directly  at  a  fully  illuminated  moon.  The  brightest  star  with  a 
wavelength  larger  than  0.5  in  has  a  spectral  irradiance  of  about 
5  X  10  ®  W/m^/m.  Inserting  this  value  into  (4.17)  and  (4.20)  results 
in 

<N>„  *  6.16  X  10  ^/M  counts/time  slot/star  , 

D 


a  value  small  enough  to  neglect  even  from  many  stars  of  the  same 
spectral  irradlance. 


Background  counts  can  also  occur  by  dark  counts  emitted  from  the 
detector.  However,  dark  counts  are  usually  on  the  order  of  10  to  20 
counts/sec.  With  the  high  data  rates  under  consideration,  the  dark 
counts  per  time  slot  are  extremely  small  and  also  negligible. 

The  above  calculations  show  that  the  worst  possible  communication 
situation  is  encountered  by  a  satellite  in  geosynchronous  orbit  pointed 
directly  towards  a  cloud-covered  section  of  the  earth  that  is  illumi¬ 
nated  by  the  sun.  Using  1  to  3  background  counts/ time  slot  for  the 
nominal  case  and  7  to  10  background  counts/ time  slot  for  the  worst  case 
would  effectively  represent  the  range  of  expected  operating  conditions. 
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CHAPTER  5 

SIMULATION  OF  THE  RECEIVER 


The  computer  simulation  of  the  receiver  is  comprised  of  two  main 
parts.  The  first  is  the  modeling  of  the  photodetector  which  generates 
the  arrival  time  records,  and  the  second  is  the  simulation  of  the  PLL. 

To  generate  the  arrival  time  records  we  follow  the  procedure 
described  in  [11].  We  know  that  the  interarrival  times  of  a  homogeneous 
Poisson  process  with  intensity,  X^,  are  independent  and  identically 
exponentially  distributed  also  with  parameter  Xp.  Using  this  property 
we  first  generate  a  sequence  of  independent  and  uniformly  distributed 
random  variables  u^ ,  U2»  ...»  Uj^  on  the  interval  (0,1).  Then  using  the 
transformation 

■  -lnCu^)  (5.1) 


a  sequence  of  Independent,  exponentially  distributed  random  variables  of 
unit  parameter  is  obtained.  To  generate  a  homogeneous  Poisson  process 
we  assign  as  the  arrival  time  of  the  first  photon,  and  then  define  ^ 
as  the  Interarrival  time  of  the  second  photon  and  as  the  interarrival 
time  of  the  kth  photon.  Now,  the  arrival  time  of  the  kth  photon,  x^,  is 


*k 


(5.2) 
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Next  a  random  PPM  intensity  sequence ,  X^t),  is  created.  The  final 
arrival  time  sequence  used  in  (2.5)  is  formed  by  solving 


*k 


/  A^z)  dz 


(5.3) 


for  the  t^'s  (11).  Appendix  A.l  contains  a  listing  of  the  Fortran  pro¬ 
gram  ARRIVAL  used  to  generate  the  arrival  time  record  for  an  tt-ary  PPM 
signal  using  the  described  method. 

The  PLL  was  simulated  using  the  Advanced  Continuous  Simulation 
Language  (ACSL).  This  language  is  designed  for  modeling  and  evaluating 
the  performance  of  continuous  systems  that  are  described  by  time 
dependent,  nonlinear  differential  equations  and  transfer  functions. 
ACSL  allows  direct  transformation  of  the  mathematical  diagram  of 
Figure  3.2  into  a  suitable  simulation  model  program.  A  listing  of  the 
ACSL  simulation  program  for  the  PLL  model  is  in  Appendix  A. 2.  This 
program  takes  an  externally  prepared  arrival  time  sequence,  processes  it 
through  the  filter  and  squaring  loop,  corn*  cts  the  shot  noise,  applies 
the  resulting  signal  to  Che  PLL,  and  gathers  data  on  the  phase  error. 
The  simplicity  of  ACSL  is  evident  in  lines  30-35  of  the  simulation 
program.  These  six  lines  contain  all  the  model  information  necessary  to 
implement  the  PLL.  In  fact,  the  filter  is  completely  described  on  line 
33  by  the  ACSL  statement  REALPL  which  implements  a  simple  low-pass 
filter.  The  integrating  effect  of  the  VCO  is  contained  in  line  35  by 
use  of  the  INTEG  statement  which  performs  an  integration  on  the  VCO 
input  control  voltage.  Most  of  the  remaining  program  is  devoted  to 
simulation  of  the  preprocessing  components  and  shot  noise  construction. 


To  show  the  validity  of  the  simulation  model,  we  present  results  of 
two  simulations:  one  without  signal  preprocessing  and  the  second  with 
preprocessing.  From  the  results  in  Chapter  2,  we  expect  phase-lock 
behavior  only  for  the  preprocessed  case  due  to  the  presence  of  the 
dominating  spectral  frequency  component  in  the  PLL  input.  Figures  5.1 
and  5.2  show  the  Fourier  transforms  of  the  PLL  input  and  VCO  output 
signal  when  preprocessing  is  by-passed,  and  where  the  frequency  scale  is 
normalized  to  the  time  slot  rate.  Note  the  lack  of  a  dominating 
frequency  in  the  input  spectrum  and  that  the  VCO  output  spectral 
component  occurs  below  the  slot  rate  at  u>^,  the  VCO  free-running 
frequency.  Figure  5.3  is  a  plot  of  Aut  versus  time.  Recall  that 
Aut**Kovc(t)t  is  the  phase  change  of  the  VCO  due  to  the  control  voltage 
applied  to  its  input.  This  change,  when  added  with  the  phase  due  to  the 
free-running  VCO  frequency,  should  equal  the  phase  of  the  input  signal. 
The  frequency  difference  between  the  slot  rate  and  u»f  was  0.103  rad/sec, 
which  should  equal  the  slope  of  Figure  5.3.  Instead,  Figure  5.3  shows  a 
very  random  behavior  rather  than  a  well-defined  slope  of  0.1  rad/sec. 
Clearly  with  no  signal  preprocessing,  and  thus  no  dominating  spectral 
component,  the  loop  cannot  achieve  phase  lock. 

Figures  5. A  and  5.5  display  the  Fourier  transforms  of  the  input  and 
VCO  output  when  preprocessing  is  used.  Again,  the  frequency  scales  are 
normalized  to  the  slot  rate.  Notice  the  very  strong  slot  frequency 
component  in  the  PLL  input  signal.  Now  the  VCO  output  spectrum  in 
Figure  5.5  has  shifted  from  its  free-running  frequency  to  match  the  PLL 
input  signal  frequency  component  at  the  slot  rate.  The  plot  of  Aut 
versus  time  of  Figure  5.6  now  has  the  desired  slope  of  0.1  rad/sec. 
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NORMALIZED  FREQUENCY  (a>/ous) 

Phase-locked  loop  Input  signal  frequency  spectrum  when  no 
preprocessing  Is  used.  Note  the  lack  of  a  dominating  fre 
quency  component. 
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Figure  5.4.  Phase-locked  loop  Input  signal  frequency  spectrum  when 

preprocessing  is  used.  The  dominating  spectral  component 
occurs  at  the  slot  frequency. 
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Figure  5.5.  VCO  output  spectrum  when  preprocessing  is  used.  Note  chat 
the  spectral  component  has  shifted  up  to  the  slot  frequency 
from  its  free  running  frequency. 


SECONDS 


do jt  versus  time  when  preprocessing  Is  used.  Note  the  slope  of 
approximately  0.1  rad/sec  Indicates  tracking  performance  of 
the  phase-locked  loop. 


CHAPTER  6 


COMPUTER  SIMULATION  RESULTS 


Computer  simulation  of  the  timing  subsystem  was  conducted  under 
varying  signal  and  background  strengths  to  represent  a  wide  range  of 
operating  conditions.  In  this  chapter  we  present  the  simulation  results 
and  compare  them  to  a  current  theoretical  expression. 

6.1  Simulation  Performance  Parameters 

The  parameter  of  Interest  is  the  phase  error,  9(t),  defined  as 


9<t)  -  0s(t)  -  0Q(t) 


(6.1) 


where,  again,  0g(t)  and  0Q(t)  are  the  PLL  input  and  VCO  output  signal 

phases,  respectively.  Once  the  loop  is  in  lock,  we  expect  the  phase 

error  to  vary  about  some  constant  value.  Assuming  9(t)  is  approximately 

Gaussian,  the  best  system  performance  indicator  of  the  PLL  is  the  RMS 

phase  error,  o,(rad). 

9 

Derivation  of  an  exact  expression  for  o  is  lengthy  and  difficult. 

9 

However,  C.  Chen  has  formulated  the  following  approximation  of  a,  for 

9 

the  receiver  timing  subsystem  presented  above.  From  [12], 


\t^\  \h[m  iV«> 

oo 

x-i-/  p(  o»* )  P*(«")  P*(u*- 
4ir  — • 


(6.2) 


xU-  /  p(o»')  p*(w")  P*(<u'-u")  H(u')  H(urw')  H*( uru")  H*(o»")  dM'dto' 


where 

K  -  expected  signal  count 
s 

H^(  to)  «  closed-loop  frequency  response  of  the  PLL 

P(oj)  *  Fourier  transform  of  received  pulse  shape 

H(o>)  -  frequency  response  of  the  detector  and  linear  filter 

H  -  number  of  time  slots/word 

T  -  slot  width, 

s 

An  important  feature  of  (6.2)  is  that  a  is  inversely  proportional  to 

9 

the  square  root  of  the  average  signal  count,  K  .  This  Important 

s 

property  provides  a  convenient  means  of  comparing  the  theory  to  the 
actual  simulation  behavior  in  the  next  section. 


6.2  Numerical  Results 

A  total  of  26  computer  simulations  were  completed  under  varying 
signal  and  background  count  conditions  for  both  4  slot  and  8  slot  PPM 
signals.  The  following  is  a  summary  of  system  parameters  used  in  the 


trials: 


time  slot  width 
time  slot  rate 
free  running  frequency 
phase  detector  gain 


T  -  1  sec 
s 

o>  -  2ir  rad/sec 
8 

<uf  -  6.18  rad/sec 
K.  -  0.3333  V/rad 


VCO  gain  constant 
loop  filter  bandwidth 


K  “1  rad/sec/V 
o 

1/t  ■  0.3333  rad/sec 
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Each  trial  represents  an  8000  sec  operation  time,  corresponding  to  2000 
data  words  (4000  bits)  for  4  slot  PPM  and  1000  data  words  (3000  bits) 
for  8  slot  PPM.  The  phase  error  was  calculated  at  the  beginning  of  each 
time  slot,  and  these  values  were  used  to  calculate  the  RMS  phase  error 
over  the  simulation  time.  Plots  of  o^  versus  time  for  each  trial  appear 
in  Appendix  B.  Notice  as  the  average  background  count,  increases 
and  the  average  signal  count,  ,  decreases,  the  RMS  phase  error 
increases  and  exhibits  progressively  erratic  behavior. 


The  a.  versus  time  plots  of  Appendix  B  are  summarized  in  Table  6.1. 
? 

Table  6.1  lists  K^,  Kg,  o^,  as  well  as  the  average  signal  count/bit, 

Kblt,  so  4  and  8  slot  PPM  formats  can  be  correctly  compared  at  the  bit 

level.  On  a  signal  count/bit  comparislon,  8  slot  PPM  has  a  higher  o, 

P 

(and  thus  inferior  performance)  than  4  slot  PPM  with  the  same  K^it  and 
K^.  Figure  6.1  shows  that  the  difference  in  performance  increases  as 
Kbit  increases.  Also,  two  of  the  simulation  trials  for  8  slot  PPM 
(Ks«*45,  Kb“5,  and  Kg«90,  Kb“20)  exhibited  cycle-slipping  behavior, 

indicating  that  8  slot  PPM  may  be  more  sensitive  to  input  signal 


variations. 


Figure  6.2  is  a  log-log  plot  of  o  versus  average  signal  count  for  4 

P 

slot  PPM  and  a  constant  of  1,  5,  and  10,  respectively.  Figure  6.3  is 

a  similar  plot  for  8  slot  PPM  with  only.  As  indicated  in 

Section  6.1,  we  expect  o.  to  be  Inversely  proportional  to  the  square 

♦ 

root  of  K  ,  corresponding  to  an  ideal  slope  of  -0.5  in  Figures  6.2  and 

3 

Each  of  the  lines  drawn  in  Figures  6.2  and  6.3  represents  a  least 
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TABLE  6.1 


NUMERICAL  RESULTS  OF  SIMULATION  TRIALS 


Average 

Signal 

Count 

K 


Average 

PPM 

Average 

RMS 

Signal 

Signaling 

Background 

Phase 

Count/Bit 

Order 

Count 

Error 

H* 

ft 

M 

*b 

<3^  (rad) 

15 

4 

1 

0.2926 

15 

4 

5 

0.3399 

20 

4 

1 

0.2609 

20 

4 

5 

0.2737 

20 
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10 

0.3146 
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1 
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Figure  6.1.  RMS  phase  error  versus  average  signal  count/bit  for  M 
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squares  fit  of  data  points  with  respect  to  a  constant  value  of  K^.  In 
Figure  6.2  the  slopes  for  4  slot  PPM  trials  are:  -0.40  for  K^l,  -0.60 
for  and  -0.58  for  K^»10.  For  8  slot  PPM  trials  in  Figure  6.3,  the 

slope  was  -0.425  for  K^-5.  These  results  compare  favorably  to  the  ideal 
slope  of  -0.5  predicted  by  (6.2).  Currently,  a  factor  of  between  three 
and  four  separates  the  theoretical  and  simulation  values.  It  is 
reemphasized  that  (6.2)  is  an  approximation,  not  an  exact  expresssion. 

The  exact  expression  for  o.  is  unwieldy  and  difficult  to  evaluate. 

9 

Efforts  are  continuing  to  obtain  a  more  accurate  approximation  of  the 
exact  expression  for  a.  [12]. 

Figure  6.4  shows  the  relationship,  for  4  slot  PPM,  between  <j^  and 
when  Kg  is  held  constant  to  40,  50,  60,  and  80  average  signal  counts. 
Figure  6.5  is  a  similar  plot  for  8  slot  PPM  for  Kg-90  only.  In  all 
cases,  the  relationship  between  the  RMS  phase  error  and  background  count 
is  not  a  simple  one.  For  lower  background  counts,  the  RMS  phase  error 
rises  moderately  with  increasing  K^.  As  the  background  count  continues 
to  Increase,  the  RMS  phase  error  rises  sharply  at  about  K^-5,  indicating 
a  rapid  deterioration  of  timing  performance. 

The  ratio  between  the  timing  offset,  AT,  and  the  time  slot  width, 
Tg,  called  the  time  slot  error  (TSE),  is  related  to  the  RMS  phase  error 
by 


tse  «  Jp-  «  . 

Ts  2* 


(6.3) 


gained  in  block  coding  (when  the  probability  of  word  error  (PWE)  is 
converted  to  PBE  at  normalized  rates)  are  almost  eliminated  as  the 
timing  error  (or  TSE)  increases  (4],  As  a  typical  example,  consider  the 
4  slot  simulation  run  of  Ks«60  and  Kb»l  of  Table  6.1.  The  RMS  phase 
error  is  0.2235  rad/sec,  corresponding  to  a  TSE  of  0.0356.  This 
indicates  the  receiver  timing  is  off  by  3.562.  The  8  slot  simulation 
trial  of  Kg-90  and  Kb*l  has  a  RMS  phase  error  of  0.2978  rad/sec  and  a 
TSE  of  0.0474,  corresponding  to  a  4.742  timing  error.  Detailed 
discussion  and  analysis  of  timing  error  effects  on  optical  PPM 
communication  systems  can  be  found  in  [1]  and  [4]. 

This  concludes  the  numerical  results  of  the  computer  simulation 


CHAPTER  7 


CONCLUSION 


In  this  report  we  have  Investigated  the  timing  performance  of  a 
phase-locked  loop  (PLL)  In  an  optical  pulse  position  modulation  (PPM) 
communication  system.  It  was  shown  that  the  output  from  a  photodetector 
alone  contains  no  dominating  spectral  frequency  component  and  is 
insufficient  as  a  signal  source  for  the  PLL.  A  simple  filter  and 
squaring  loop  was  shown  to  generate  a  strong  spectral  component  at  the 
slot  frequency.  Preprocessing  the  signal  in  this  manner  enables  the  PLL 
to  lock  up,  track  the  signal,  and  provide  timing  information  to  the 
receiver  decision  section. 

We  have  discussed  various  PLL  design  considerations.  Once  design 
parameters  are  fixed,  a  change  in  the  input  signal  magnitude  can  alter 
the  PLL  performance.  To  alleviate  this  problem,  an  expression  for  an 
attenuation  constant  was  derived  which  would  limit  the  average  PLL  input 
signal  to  unity,  allowing  the  intended  design  parameters  to  remain 
unchanged  over  a  wide  range  of  signal  inputs. 

A  general  satelllte-to-satellite  link  equation  was  derived  to 
calculate  the  average  number  of  signal  counts  emitted  by  the  receiver 
detector.  Average  signal  count  depends  on  many  factors,  including 
physical  component  values  (aperture  size,  wavelength,  etc.),  data  rate, 
and  pointing  error.  Sources  of  background  radiation  such  as  the  moon. 


stars,  and  earth  were  examined  and  It  was  shown  that  a  cloudy. 
Illuminated  section  of  the  earth  results  In  the  largest  amount  of 
background  radiation  at  the  receiver.  Typical  component  values  for  the 
link  were  chosen,  and  average  signal  and  background  counts  calculated. 

Finally,  the  entire  timing  subsystem  was  simulated  on  a  digital 
computer  and  the  RMS  phase  error  of  the  PLL  was  calculated.  The 
simulations  showed  the  RMS  phase  error  Is  (approximately)  Inversely 
proportional  to  the  square  root  of  the  average  signal  count  for  both  4 
slot  and  8  slot  PPM  signals.  This  agrees  with  a  theoretical  expression 
recently  derived  for  this  timing  subsystem  which  predicts  an  exact 
inverse  relationship  between  the  RMS  phase  error  and  the  square  root  of 
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PROGRAM  LISTINGS 

A.l  Arrival 

i 

1 

PROGRAM  ARRIVAL 

1 

2  C 

t 

3  C 

THIS  PROGRAM  CONSTRUCTS  THE  PHOTON  ARRIVAL 

4  C 

TIMES  FOR  A  M-ARY  PPM  DATA  SEQUENCE  DETECTED 

5  C 

BY  A  PHOTODETECTOR. 

-  . 

6  C 

PULSHI-AMPLITUDE  OF  SIGNAL  PULSE 

ft 

7  C 

NSLOT-ORDER  OF  THE  PPM  SIGNAL 

1 

8  C 

TMAX-LENGTH  OF  RECORD  DESIRED  IN  SECONDS 

1 

9  C 

DARK-BACKROUND  RADIATION  PER  PULSE 

10  C 

TSLOT-DURATION  OF  A  SINGLE  PULSE 

11  C 

12 

DIMENSION  ARIVAL( 1024) 

13 

INTEGER  ISLOT(25) 

l 

14 

DATA  NSLOT/4/ .TMAX/9000.0/ 

4m — — 

I 

15 

DATA  PULSHI/40./, DARK/10. /.TSLOT/l.O/ 

V 

- 

16 

AREA  -  0.0 

17 

TNOW  -  0.0 

18 

TIMLST  -  0.0 

19 

CALL  RANSETC31) 

: 

20 

DO  10  1-1,10 

— 

— . 

1 

21  10 

U  -  RANFQ 

1 

22 

CALL  DATGEN(  I SLOT ,  NSLOT) 

23 

KK  -  1 

24  100 

CALL  POISSN(ARIVAL, TIMLST, 1024, 1.0) 

H 

25 

DO  200  I  -  1,1024 

B 

ft 

26  170 

DIFF1  -  (ISLOT(KK)  *  PULSHI  +  DARK)  *  TSLOT 

_ 

■ 

1 

27 

DIFF2  -  ARIVAL(I)  -  AREA 

28 

IF  (DIFF2  .LE.  DIFF1)  GO  TO  180 

29 

KK  -  KK+1 

30 

TNOW  -  TNOW  +  TSLOT 

31 

AREA  -  AREA  +  DIFF1 

32 

IF  (KK.LE. NSLOT)  CO  TO  170 

i 

33 

CALL  DATGEN(ISLOT, NSLOT) 

34 

KK  -  1 

35 

GO  TO  170 

36  180 

ARIVAL(I)  -  TNOW  +  DIFF2/DIFF1 

. 

37  200 

CONTINUE 

38 

WRITE(3)  (ARIVAL(IW),IW- 1,1024) 

i 

39 

IF  (TNOW  .LT.  TMAX)  GO  TO  100 

40 

RETURN 

L 

41 

END 

• 

42  C 

43  C 

l 

i 

— 

r 

1 
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44  C 

45  SUBROUTINE  POISSN( ARIVAL , TIMLST , NUMTRY , PARAM) 

46  C 

47  C  SUBROUTINE  POISSN  GENERATES  THE  POISSON  INPUT 

48  C  SEQUENCE  OF  UNIT  PARAMETER  USING  UNIFORM 

49  C  RANDOM  VARIABLES. 

50  C 

51  DIMENSION  ARIVAL(l) 

52  1  UNIF  -  RANF() 

53  IF  (UNIF. EQ. 0.0)  GO  TO  1 

54  ARIVAL(l)  -  TIMLST  -  (AL0G( UNIF) /PARAM) 

55  DO  10  J- 2, NUMTRY 

56  12  UNIF  -  RANF() 

57  IF  (UNIF. EQ. 0.0)  GO  TO  12 

58  ARIVAL(J)  -  ARIVAL(J-l)  -  (ALOG( UNIF) /PARAM) 

59  10  CONTINUE 

60  TIMLST  -  ARIVAL( NUMTRY) 

61  RETURN 

62  END 

63  C 

64  C 

65  SUBROUTINE  DATGEN( ISLOT, NSLOT) 

66  C 

67  C  SUBROUTINE  DATGEN  GENERATES  THE  RANDOM 

68  C  SIMULATED  PPM  DATA  SEQUENCE. 

69  C 

70  INTEGER  ISLOT(l) 

71  INTEGER  MASK, IRAND 

72  MASK  -  NSLOT 

73  DO  10  1-1, NSLOT 

74  10  ISLOT(I)  -  0 

75  IRAND  -  IFIX(RANF()  *  1000) 

76  IRAND  -  MOD( IRAND, MASK) 

77  ISLOT ( IRAND+1 )  -  1 

78  RETURN 

79  END 


A. 2  Simulation  of  a  Phased-Locked  Loop 


1  PROGRAM  SIMULATION  OF  A  PHASE-LOCKED  LOOP 

2  “THIS  PROGRAM  IS  DONE  IN  THE  ADVANCED" 

3  "CONTINUOUS  SIMULATION  LANGUAGE  (ACSL)“ 

4  •  * 

5  INTEGER  INDEX 

6  REAL  W0,KV 

7  REAL  ARIVAL( 2049) 

8  CONSTANT  THETA  -  0..KV  -.3333,... 

9  FIC2  -  0.0, 

10  PHIO  -  0.0, WO  -  6.18,  ... 

11  INITIAL 

12  CINTERVAL  CINT  -  0.125 

13  NSTEPS  NSTEP  -  3 

14  IALG  -  4 

15  INDEX  -  1 

16  ENDTIM  -  8050.0 

17  PROCEDURAL  (ARIVAL-ARIVAL) 

18  CALL  SETVAL(ARIVAL) 

19  END  $"OF  PROCEDURAL" 

20  END  $  "  OF  INITIAL  " 

21  "  " 

22  DYNAMIC 

23  DERIVATIVE  $“THIS  SECTION  CONTAINS  THE  PLL" 

24  "SIMULATION  MODEL  EQUATIONS.” 

25  • 

26  "  " 

27  PROCEDURAL  (SHOT-T.ARIVAL, INDEX) 

28  CALL  SHOTNZ(T,ARIVAL, INDEX, SHOT) 

29  ENDS  "  OF  PROCEDURE  " 

30  PHASE  -  WO*T+PHI+THETA 

31  COPHAS-COS( PHASE) 

32  MIXOUT  -  SHOT*COS( PHASE) 

33  FILOUT  -  REALPLC 3., MIXOUT, 0.0) 

34  PHIDOT  -  KV*FILOUT 

35  PHI  -  INTEG( PHIDOT, PHIO) 

36  TERMT(T.GE. ENDTIM) 

37  END  $  "  OF  DERIVATIVE  ” 

38  WRITE(l)  PHASE  $  "SAVES  THE  VCO  OUTPUT  PHASE  FOR  ANALYSIS' 

39  END  $  "  OF  DYNAMIC  SECTION  " 

40  END  $  "  OF  PROGRAM  " 
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43 

44  C 

45  C 

46  C 

47  C 

48  C 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58  C 

59  C 

60  C 

61  C 

62  C 

63  C 

64 

65 

66 

67  10 

68 

69 

70 

71  20 

72 

73 

74 

75  100 

76 

77 

78  101 

79 

80 

81  C 

82  C 

83  C 

84  C 

85  C 


SUBROUTINE  SETVAL(ARIVAL) 

SETVAL  READS  IN  THE  FIRST  TWO  BLOCKS  (2048  VALUES)  OF 
THE  ARRIVAL  TIMES  FROM  THE  EXTERNALLY  PREPARED  ARRIVAL 
TAPE. 

REAL  ARIVAL(2049) 

READ(3)  ( ARIVAL( J) ,  J-l , 1024) 

READ(3)  (ARIVAL(J),J-1025,2048) 

RETURN 

END 


SUBROUTINE  SHOTNZ(T,ARIVAL, INDEX, SHOT) 

SUBROUTINE  SHOTNZ  EVALUATES  THE  SHOTNZ  SIGNAL 
FOR  INPUT  TO  THE  PLL  AT  A  GIVEN  TIME  T.  SHOTNZ 
RETURNS  THE  CURRENT  VALUE  IN  THE  PARAMETER  SHOT. 


REAL  ARIVAL( 2048 ) 

INTEGER  INDEX 
Q-T+10 

IF  (ARIVAL( INDEX). GE.Q)  GO  TO  100 
INDEX  -  INDEX  +  1 
IF  (INDEX.LT. 2049)  GO  TO  10 
DO  20  1-1,1024 

ARIVAL(I)  -  ARIVAL( 1+1024) 

READ(3)  (ARIVAL( I) ,1-1025,2048) 

INDEX  -  1025 
GO  TO  10 
SHOT  -0.0 
FLAG  -  1 
II  -  0 
II  -  II  +1 

SHOT  -  SHOT  +  GATE( Q-ARIVAL( INDEX-II) , IFLAG) 
IF  (IFLAG.GE.O)  GO  TO  101 

THE  NEXT  LINE  TAKES  THE  SHOT  NOISE  VALUE  AND 
SQUARES  IT  AND  DIVIDES  BY  THE  ATTENUATION 
CONSTANT  OF  THE  PREPROCESSING  SECTION. 
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86 

87 

88 

89 

90  C 

91  C 

92 

93  C 

94  C 

95  C 

96  C 

97 

98 

99 
100 
101 
102 

103 

104  10 

105 

106 
107 


SH0T1  -  SH0T*SH0T/ 40.528 
SHOT  -  SH0T1 

RETURN 

END 


REAL  FUNCTION  GATE(X, IFLAG) 

FUNCTION  GATE  SIMULATES  THE  COMBINED 
PREPROCESSING  FILTER  IMPULSE  RESPONSE. 

REAL  PI2 

DATA  PI2/6. 283185307/ 

IF  (X.GT.1.0)  GO  TO  10 
G  -  SIN(PI2*X) 

IFLAG  -  1 
GATE  -  G 
RETURN 

IFLAG  -  -1 
GATE  -  0.0 
RETURN 
END 


RMS  PHASE  ERROR  (RAD) 
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RMS  PHASE  ERROR  VERSUS  TIME  FOR  SIMULATION  TRIALS 


Figure  B.l.  Simulation  trial  for  Ks  -  30,  Kj,  «  1  and  M  -  4. 
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Figure  B.3.  Simulation  trial  for  K  »  40,  K.  -  1  and  M 
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SECONDS 


Figure  B.5.  Simulation  trial  for  Kg  -  40,  ■  10  and  M  ■  4. 
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Figure  B.7.  Simulation  trial  for  K  -  50,  K.  «  5  and  M  ■ 
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Figure  B.8.  Simulation  trial  for  K  -  50,  K.  -  10  and  M  ■  4 
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Figure  B.10.  Simulation  trial  for  K  -  60,  K. 
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2000.  4000.  6000. 

5000.  5000.  7000. 

SECONDS 


8000 


Figure  B.ll.  Simulation  trial  for  K  *60, 


*  5  and  M  *  4 
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Figure  B.12.  Simulation  trial  for  Kg  »  60,  *10  and  M  ■  4. 
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SECONDS 

Figure  B. 16.  Simulation  trial  for  Kg  •  80,  «  10  and  M  *  4. 
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Figure  B.18. 
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Figure  B.19.  Simulation  trial  for  K  -  75 
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Figure  B.22.  Simulation  trial  for  Kg  -  90,  -  1  and  M  -  8. 
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Figure  B.24.  Simulation  trial  for  Kg  -  90,  ■  10  and  M  ■  8. 
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